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Abstract. We develop further a kinetic theory of strings and derive a transport equation 
for a network of cosmic strings with Nambu-Goto evolution, interactions and background 
gravitational effects taken into account. We prove an H-theorem and obtain necessary and 
sufficient conditions for a thermodynamic equilibrium. At the lowest order the equilibrium is 
estimated by the von Mises-Fisher distributions parametrized by mean directions and disper- 
sions of the right- and left-moving tangent vectors. Under assumption of a local equilibrium 
we derive a complete set of hydrodynamic equations that govern the evolution of strings 
on large scales. We also argue that on small scales the assumption of a local equilibrium 
would break down, and non-equilibrium steady states, described by the Sinai-Ruelle-Bowen 
distributions, should be used instead. 
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1 Introduction 

Cosmic strings are predicted by many models of symmetry breaking phase transitions [1] and 
give rise to very distinct and detectable signatures such as gravitational lensing [2], CMB 
non-Gaussianities [3], gravitational waves [4], ultrahigh energy cosmic rays [5], radio signals 
[6] etc. It is also believed that cosmic superstrings can be formed at the end of brane inflation 
[7] which opens a possibility of testing the models of string theory in the cosmological settings. 
However, all of the phenomenological predictions are based on the statistical properties of 
cosmic strings, which remain largely unexplored. 

In the early work on strings statistics [8], it was expected that a typical length of closed 
loops would scale linearly with time, but some old numerical simulations showed no evidence 
for such scaling [9-13]. The average distance between long strings was observed to scale 
linearly with time, but the size of the smallest wiggles as well as the size of loops remained 
at the resolution scale of the simulations. More recently, the high-precision numerical simu- 
lations [14] with much finer resolution were developed to demonstrate that the long strings 
do obey a scaling law on a wide range of scales described by a universal power spectrum. 
Moreover, there is evidence suggesting that after a long transient regime, characterized by 
an excessive production of small loops at the scale of initial conditions, the characteristic size 
of large loops starts to grow linearly with time [15]. (See also [16] for the numerical results 
and interpretations reported by other groups.) Unfortunately, the dynamical range of these 
simulations is not sufficiently large to confirm the numerical results with a high confidence 
level. 

In addition to numerical studies, a number of analytical models were proposed [17, 18], 
but the downside of most models is that they contain phenomenological parameters, such 
as the mean velocity of strings. To overcome these problems and to better understand 
various nonlinear stringy phenomena, such as cross correlations [19] and semi-scaling [20], we 
developed a kinetic theory of strings [21]. The stochastic dynamics of strings was modeled 
with a transport equation which is capable of describing long strings with Nambu-Goto 
evolution as well as interactions taken into account; and the corresponding system of coupled 
transport equations which can simultaneously describe long strings as well as closed loops. 
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In this paper we extend the kinetic theory approach further and derive a fully hydrodynamic 
description of strings. 

The paper is organized as follows. In the next section we derive a simplified version of the 
transport equation for strings. In Sec. 3 we prove an i?-theorem and derive an equilibrium 
condition for strings. In Sec. 4 we prove a conservation theorem and obtain a general set 
of conservation equations. The hydrodynamic equations for von Mises-Fisher distributions 
are derived in Sec. 5 under assumption of a local thermodynamic equilibrium. In Sec. 6 
we relax the equilibrium assumption and discuss the consequences of nonequilibrium steady 
states. The main results of the paper are summarized in Sec. 7. 

2 Transport Equation 

The dynamics of cosmic strings without interactions, such as reconnections of nearby strings 
or production of string loops, is described fairly well by Nambu-Goto action, 

Sng = y V-^et{Kb)dadt (2.1) 

where ^ is the string tension and 

hab = gf^uX'^aX'^h- (2.2) 

The corresponding equation of motion is 

x(^,^'^ + ri;>-V,<, = o. (2.3) 

After gauge fixing (/iqi = 0, /iqo + hu = 0, t = x^) on a flat background we get a system of 
equations for a three-vector, x, 

(fx (fx 
ax \ / \ 



4tJ ' V*?j 

(ix dx ^ ^2 

da dt ' 

where the speed of light is set to one. The solution of (2.4) can be described in terms of 
right-moving a(o" — t) and left-moving h{a + t) waves 

x(a,t) = "("-^^+^(" + ^^ (2.7) 

with additional constraints, 

1. (2.8) 



da 




dh 


da 




da 



due to (2.5) and (2.6). 

In general the exact evolution of a given infinitesimal segment of string is known only if 
all of the functions, Si{a) and h(a), for all strings are known. In the kinetic theory of strings 
such evolution is described by a system of equations analogous to the Bogoliubov-Born- 
Green-Kirkwood-Yvon hierarchy for particles [21]. Truncating the hierarchy at a given order 
is equivalent to neglecting higher-order effects and is often a starting point in a derivation of 
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a transport equation. The first attempt to drive a transport equation for strings was given 
in Ref. [21] where the state of the system was modeled with a 9 + 1 dimensional distribution 
function /(A, B, x, t) of the correlation vectors A and B, comoving position x, and conformal 
time t. In this paper we make a crucial simplification by concentrating instead on a 7 + 1 
dimensional distribution /(A,B,x, t) where 

J da , ~ db , , 

A = and B = — 2.9 

da da 

are the unit three-vectors corresponding to tangent vectors of right- and left-moving waves. 
In a Priedmann universe described by a Friedmann-Robertson- Walker metric, 

ds'^ = a{t)^idf -dx'^), (2.10) 

one can show that the unit vectors evolve according to 

|a = -?^(B-(A.B)A) (2.11) 

^B = -?^(A - (A • B)B) (2.12) 

and the energy density decays as, 

^; = -«1(A + b)V. (2.13) 

where H = ^ is the Hubble rate in conformal coordinates (See, for example, Ref. [11]). Then 
one can derive a transport equation for /(A, B, x, t) by following the evolution of /(A, B, x, t) 
from time t to time t + 6t. For an infinitesimally small 5t, we can use (2.11), (2.12) and (2.13) 
to express the distribution at time t + 5t in terms of the distribution at time t, 

f (^A - n (b - (^A ■ St,B - n (^A - (^A ■ Bj Bj Jt, x + ^ (^A + b) 5t, t + 6t 



l-n-[A + B) j/(A,B,x,t)+ (^^j Jt. (2.14) 

The "collision" term on the right-hand side represents the Nambu-Goto evolution and inter- 
actions. By expanding to the linear order in 6t, we obtain the following equation 

where the dependence on the comoving position x and conformal time t is suppressed for 
brevity of notations. Note, that the use of the simplified distribution function /(A,B,x, i) 
does not necessarily imply that the correlations must be small. In fact, the long-distance 
correlations may still exist whenever the distributions of A and B are peaked sharply around 
a particular direction. 

The effects of "collisions" can be estimated as in Ref. [21]. The Nambu-Goto "collisions" 
proceed with the same rate, regardless of the orientations of string world sheets, as the right- 
and left-moving waves pass through each other. The rate of such "collisions" is given by 

revolution ~ 1 , (2.16) 
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where ^min describes the resolution scale of the kinetic theory. Since the Nambu-Goto action is 
only accurate on sufficiently large scales there is a lower bond Zmin ^ fi~^^^. However, on the 
cosmological settings the resolution scale might also be set by the scale of the initial conditions 
or by the scale of the gravitation backreaction. The exact value of /min is unimportant for 
the subsequent discussion. 

The rate of the second kind of "collision", due to interactions between a segment de- 
scribed by tangent vectors A and B with another segment described by tangent vectors A' 
and B', is given by a probability of intersecting in the target space. In a universe with three 
large spatial dimensions the probability must be proportional to the local energy density, p, 
the intercommutation probability, p,^ and also to what we call a cross-sectional volume. Intu- 
itively, the cross-sectional volume describes a likelihood of interactions of string segments due 
to relative orientations of their world sheets (1, A) A (1, B) and (1, A') A (1, B'). Altogether 
we get the following expression for the rate of interactions, 



1, A) A (1, B) A (1, A') A (1, B') (2.17) 
dAdB/(A,B) (2.18) 



Tintcraction ~ TP P 

where 



and 7 is some constant. 

When a "collision" takes place the corresponding opposite-moving tangent vectors are 
interchanged, e.g. 

(A',B) and (A,B') (A,B) and (A',B'), (2.19) 

or 

(A,B) and (A',B') ^ (A',B) and (A,B'). (2.20) 

This gives rise to a couple of integral terms which can be combined together to yield the 
following expression, 

df\ f.X,.^nr^ , ,/(A',B)/(A,B')-/(A,B)/(A',B') 



j dA dB (revolution ~l~ rintcraction 



C>t J coll J P 

(2.21) 

By substituting (2.21) into (2.15) we obtain the final form of the transport equation. 



^ dA'dB'7'(A,B,A',B')(/(A',B)/(A,B')-/(A,B)/(A',B')) (2.22) 
where 

P(A, B, A', B') ~ + 7 p (1, A) A (1, B) A (1, A') A (1, B') . (2.23) 

^min P 

Note that the function P(A, B, A', B') is invariant under permutations of variables. This 
property will turn out to be crucial for proving an iJ-theorem and a conservation theorem 
in the following sections. 



^The intercommutation probability is a free parameter of order one for the field theory strings but can also 
be much smaller than one of the string theory strings. 
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3 Equilibrium Distributions 

Consider a localized version of the transport equation (2.22) where the spatial inhomo- 
geneities and gravitational effects are ignored, i.e. 

^/(A, B) = J dA'dB' P(A, B, A', B') (/(A', B)/(A, B') - /(A, B)/(A', B')) . (3.1) 

If we assume that a local equilibrium (defined by ^/eq(A, B) = 0) is established much faster 
than all other time scales, then the right-hand side of (3.1) must also vanish, 

J dA'dB' ViA, B, A', B') (/(A', B)/(A, B') - /(A, B)/(A', B')) = 0. (3.2) 

We are now ready to prove an ff-theorem for strings which is analogous to the famous 
Boltzmann's i/-theorem for particles. In particular, we will show that 

H{^,t) = I dAdB/(A,B)log(/(A,B)), (3.3) 

can never increase with time 

dH , , 

^ < 0, (3.4) 

whenever /(A,B) satisfies the local transport equation (3.1). 
By differentiating (3.3) with respect to time, we obtain 



i|f = /dAdBM|a(i + ,„,(;(A,B))) (3.5) 

and after substituting (3.1) into (3.5) we arrive at the following equation 

y dAdBdA'dB' P(A,B, A',B') (/(A', B)/(A, B') - /(A, B)/(A', B')) x 



dR 
'dt 



l + log(/(A,B)) . (3.6) 



Since P(A, B, A', B') is invariant under permutations of variables, the right-hand side of 
(3.6) can be rewritten by interchanging the dummy variables of integration A, B and A', B', 

^ = J dA'dB'dAdBP(A',B',A,B) (/(A, B')/(A', B) - /(A', B')/(A, B)) x 

'l + log(/(A',B'))).(3.7) 



Note that (3.6) and (3.7) only differ in an argument of the logarithms. Thus, we can add 
(3.6) and (3.7) to get 

2^ = y" dA'dB'dAdBP(A',B',A,B) (/(A, B')/(A', B) - /(A', B')/(A, B)) x 

(2 + log(/(A,B)/(A',B')))(3.8) 



- 5 - 



and by interchanging A and A' we obtain another equivalent equation, 

2^ = y dA'dB'dAdBP(A',B',A,B) (/(A', B')/(A, B) - /(A, B')/(A', B)) x 

(2 + log(/(A',B)/(A,B'))).(3.9) 

To arrive at the final expression we add (3.8) and (3.9), 
^ = -\l dA'dB'dAdBP(A',B',A,B) (/(A, B')/(A', B) - /(A', B')/(A, B)) x 

log(/(A,B')/(A',B)) -log(/(A,B)/(A',B'))) . (3.10) 
But since log is a monotonically increasing function we conclude that 

(/(A, B0/(A', B) - /(A', B')/(A, B)) (log(/(A, B0/(A', B)) - log(/(A, B)/(A', B'))) > 

as well as ^ < 0, which is exactly the statement of the iJ-theorem for strings (3.4). Note 
that the equality ^ = is realized if and only if 

/e,(A, B')/e,(A', B) = /e,(A', B')/e,(A, B) (3.11) 

corresponding to an equilibrium distribution /eg(A,B). 

In addition, we can show that for an arbitrary equilibrium distribution, there exists a 
factorization 

feg{A,B) = p pa{A) pb(B) (3.12) 
where Pa(A) and Pb(B) are the normalized probability distributions, i.e. 

J dA pa{A) = JdB ph{B) = 1. (3.13) 

Indeed, for an arbitrary pair of unit vectors A" and B" such that /(A",B") ^ 0, we can 
define 

Pa{A)^ and p,(B) = — ^i^^^i^l^. (3.14) 

/dA'/e,(A',B") "^'^ /dB'/e,(A",B') 

Then it is easy to check that the condition (3.11) implies (3.12), 

feg{A, B')/eg(A', B) = /e,(A', B')/e,(A, B) V A, B, A', B' ^ 

rdA'dB7eo(A',B')/eo(A,B")/eg(A",B) 
p Pa(A)pb(B) = ^ . . ^.^egv . JJeqy ) ^ , ^ 

^ ' ^ ' /dA'/e,(A',B")/dB'/e,(A",BO ^ ' 

J dAMB7e,(A^ B')f,,{A",B")f,g{A, B) ^ . 

/dA'dB7e,(A',B')/e,(A",B") ^''^ ' 

Of course, the opposite is also true, and (3.12) implies (3.11). Thus, we have four equivalent 
definitions of a local equilibrium of strings 

^/e,(A, B) = 0^ B)] = 4^ /e,(A, B) = p pa{A) p,{B) ^ 

/eg(A,B')/eg(A',B) = /e,(A',B')/eg(A,B) V A,B,A',B'. (3.16) 
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4 Conservation Theorem 

Consider an arbitrary quantity Q{A, B) which is conserved during "colhsions" (Nambu-Goto 
or interactions), i.e. 

Q(A, B) + g(A', B') = Q(A', B) + Q(A, B'). (4.1) 
Then we can integrate both sides of the transport equation (2.22) to get 

/dAdB«(A,B)(| + i(A + B).| + «(i(A + B)^-A.^-B.A))y(A,B) = 

J dAdBdA'dB' Q(A, B) P(A, B, A', B') (/(A', B)/(A, B') - /(A, B)/(A', B')) , (4.2) 

where the right-hand side can be rewritten as 

j dAdBdA'dB' Q(A,B)P(A,B,A^ BO (/(A', B)/(A,B')- /(A, B)/(A',B')) = 

j dAdBdA'dB' Q(A', B) P(A, B, A', B') (/(A, B)/(A', B') - /(A', B)/(A, B')) = 

J dAdBdA'dB' Q(A, B') P(A, B, A', B') (/(A', B')/(A, B) - /(A, B')/(A', B)) = 

J dAdBdA'dB' Q(A', B') P(A, B, A', B') (/(A, B')/(A', B) - /(A', B')/(A, B)) 

due to invariance of 'P(A, B, A', B') under permutations of variables. By adding these four 
equivalent expressions together and dividing by four we obtain yet another equivalent ex- 
pression 

^ j dAdBdA'dB'P(A,B,A',B') (/(A', B)/(A, B') - /(A, B)/(A', B')) x 

(q(A, B) + Q(A', B') - Q(A', B) - Q(A, B')) (4.3) 

which must vanish for any conserved quantity Q(A,B) due to (4.1). Thus, the right-hand 
side of (4.2) is identically zero and we arrive at the following equation, 

/dldB Q(A,B) (1 + ^ (a + b) . I- + W (i (a + b)'^ - A . ^ - B . |.)) /(A.B) = 0. (4.4) 

It is now convenient to define an expectation value 

iO) = ^j dAdB 0(A,B)/(A,B). (4.5) 
and to rewrite (4.4), after integrating by parts, as 

where the summation over repeated indices is implied. In the literature on the kinetic theory 
of particles a similar expression is known as a conservation theorem. By applying the conser- 
vation theorem to different (locally) conserved quantities we can derive different conservation 
equations. 
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In particular, an equation for a local conservation of energy is obtained by setting Q = 1 
in (4.6), 

¥t^^\^/ + + (1 + A,B,) = 0. (4.7) 

Another important example is a local conservation of right-moving and left-moving tangent 
vectors. This gives us six more equations, one for each of the components of each of the 
tangent vector. By setting Q to Ai, A2, A^, Bi, B2 or ^3 in (4.6), we obtain the following 
equations 

g-^P {Ai) + -—p {AjA, + B,A,) + np {AjBjAi + Ai + B,) = 0, (4.8) 

g-^P {B,) + -—p {AjBi + BjBi) + np {AjBjB, + A, + B.,) = 0. (4.9) 

Altogether, there are seven conservation equations [i.e. (4.7), (4.8) and (4.8)], but 
this is certainly not all there is [e.g. Q(A) = 5(A — A') is conserved for an arbitrary A']. 
In fact, there is an infinite number of conservation conditions that one could derive, and 
it may seem as if the system is overconstraint. (Note that there are only five conserved 
quantities in a kinetic theory of particles: mass, energy, and momentum, which give rise 
to five hydrodynamic equations: the continuity equation, heat equation and Navier-Stokes 
equation.) However, as we will see in the following section, seven is the exactly the number 
of equations needed to describe the hydrodynamics of strings to the leading order, and all 
other conservation conditions become relevant only at higher orders in perturbation theory. 



5 Hydrodynamic Equations 

In Sec. 3 we have shown that a local equilibrium is established if and only if /(A, B) 
factors into a product of the energy density, p, and two normalized probability distribution 
functions on a unit sphere, Pa(A) and pb(B). Clearly, an arbitrary choice of the normalized 
distributions j'a(A) and Pb(B) would satisfy the equilibrium conditions (3.16) and additional 
considerations are needed to determine the most physical distributions. 

From the conservation of energy and of tangent vectors (of right- and left-moving waves) 
we derived equations (4.8), (4.9), and (4.7), which can be rewritten under assumption of a 
local equilibrium as 

WtP'^^'^ + I'&.P^^^'^^^ + ^^^^ (i?,)) + ^p((^,A,) {B,) + {A,) + {Bi)) = 0, (5.1) 

^^P{B,) + + (^*) {Aj)) + np{{BjBi) (A,) + (A,) + (Bi)) = 0, (5.2) 

|p + {{A,) + {B,)) + npil + {A,) (B,)) = 0. (5.3) 

This is a system of seven equations at the lowest order in perturbation theory, but for a self- 
consistent description, the energy distribution function /(A,B) must also be estimated with 
seven parameters. One parameter is the normalization constant, p, and the simplest three 
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parameter family of normalized distributions is given by the von Mises-Fisher distributions 
on a unit sphere [22] , 

^'"(^) = -A ■\(M \ '^^P ( I (^•^) 

47rasmh(l/a) \ a ) 



and 



^''(^^ = 4vr/3sinh(l//3) ^"^ 0^ ^'-'^ 



where u, v are the mean directions and a, (3 are the dispersion parameters of right- and 
left-moving tangent vectors. The main motivation to use the von Mises-Fisher distributions 
is that they are mathematically tractable and at the same time give a good approximation 
of the Gaussian distributions wrapped on a sphere [23] . The corresponding lowest statistical 
momenta, 

{Ai) = u,F{a), (5.6) 
(Bi) = ViF{^), (5.7) 
{AiAj) = dijaF{a) + muj (1 - 3aF{a)) , (5.8) 

and 

{B^B,) = 5i,pF{P) + UiU, (1 - 3/3F(/3)) (5.9) 
where, for convenience, we defined 

F{x) = coth -X (5.10) 

which can be approximated in the limit of small dispersions as 

F{x)^l-x. (5.11) 

To obtain a complete system of hydrodynamic equations, we can plug the expressions 
(5.6), (5.7), (5.8) and (5.9) into the equations (5.1), (5.2), (5.3). The resulting equations are 

d 1 ^ 

—puiF{a) + i^-Q^P iSijaF{a) + muj (1 - 3aF{a)) + UiVjF{a)F{P)) + (5.12) 

np {uiF{a) + ViF{l3) (1 + aF{a)) + UiUjVjF{p) (1 - 3aF(a))) = 0, 

g-^Pv^Fm + 2d^^P i^iiPPiP) + (1 - 3/3F(/3))) + v,u,F{a)Fm + (5.13) 



Up {viF{/3) + UiFia) (1 + /3F(/3)) + v^UJVJFia) (1 - 3/3F(/3))) = 0, 



and 



o 1 ^ 

+ 2 i^j^i^) + ^iFiP)) + (1 + UiVjF{a)Fm = 0. (5.14) 

There are seven equations with seven unknowns: three scalar quantities, p, a, /3, and two 
vector quantities, u, v, constrained to unity, i.e., |u| = |v[ = 1. This system of partial dif- 
ferential equations is a huge simplification in comparison to the integro-differential transport 
equations of the kinetic theory (2.22). 



-9- 



Note that the hydro dynamic equations can be simphfied further on a Minkowski back- 
ground (i.e. 7i = 0), 

d 1 ^ 

—puiF{a) + i^-Q^P {5ijaF{a) + UiUj (1 - 3aF{a)) + UiVjF{a)F{(3)) = 0, (5.15) 
-QlPViF{l3) + 2 {5i,pF{P) + v,vj (1 - 3/3F(/3))) + ViU,F{a)Fm = 0, (5.16) 

and in the hmit of small dispersions a and /3, 
d Id 

— pnj(l - q) + ("^^i" + ^'^J ~ + ^^i'yi(l - a - /3)) = 0, (5.18) 

-p^;,(l - /3) + 2d^.P^^i^^ + (1 - 3^)) + - a - /3)) = 0, (5.19) 

d 1 ^ 

di''^2d^/ ^""^'^^ " + ""^'^^ - ^)) = 0. (5.20) 

Evidently, even without gravitational effects, the evolution of right-moving waves, described 
by Ui and a, and left-moving waves, described by Vi and /3, is coupled not only through 
energy density p, but also directly though the cross terms 

^-^pviUjF{a)F{p) and ^-^puiVjF{a)F{^). (5.21) 

These terms give rise to the so-called cross correlations between opposite moving waves [19], 
which are often overlooked in the literature on cosmic strings. Our analysis shows that the 
cross correlations appear at the leading order and should not be ignored. 



6 Non-Equilibrium Steady States 

The reader might be puzzled by the fact that the hydrodynamic equations derived in the pre- 
vious section do not depend on the intercommutation probability p. At first glance this seems 
very counterintuitive, given that there are numerical simulations that have demonstrated the 
dependence of string dynamics on p, at least on small scales. Since the only assumption that 
we have made so far is the assumption of a local thermodynamic equilibrium, there must be a 
critical scale below which the assumption is no longer valid. The scale can be quite small for 
the field theory strings with p ~ 1 but can also be rather large for the string theory strings 
with p <C 1. 

Evidently, on small scales the network of strings must not behave as a closed Hamil- 
tonian system whose equilibrium distribution would be given by the Liouville measure. For 
more general hyperbolic dynamical systems the nonequilibrium steady states provide a more 
accurate description of the dynamics. These states were originally proposed by Sinai [24], 
Ruelle [25] and Bowen [26] and go by the name of the Sinai-Ruelle-Bowen (SRB) measures 
[27]. In this paper we will refer to them as the SRB distributions. Of course, an arbitrary 
local distribution can be substituted into (4.7), (4.8), and (4.9) to yield a new set of hydrody- 
namic equations for strings, and thus the main problem is to estimate the SRB distribution 
for strings. 
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To derive a local distribution of strings we will apply the so-called thermodynamic 
formalism [27] with many ideas borrowed from the conventional statistical mechanics, but 
one important difference. In the statistical mechanics we are usually interested in the phase 
space states, X, when in the dynamical system of strings the key role is played by the phase 
space trajectories, X[t). Thus, it is convenient to think of time, t, as a thermodynamical 
volume which is a conjugate variable to the topological pressure defined as 

p = lim 7^ log 2:, (6.1) 



where 



Z = j dX(0) exp i-j^Yl Xi{X{t))dt j (6.2) 



Ai>0 

is the dynamical partition function [27]. The sum of local Lyapunov exponents Xi (defined 
as a local rate of separation of the nearby trajectories) is taken over directions corresponding 
to only positive global Lyapunov exponents (defined as a rate of separation of nearby 
trajectories in the limit of infinite times). 

According to the partition function (6.2) the trajectories with large Lyapunov exponents 
are exponentially suppressed. This effect can be easily estimated for strings from an expected 
collision rate of a given configuration. In the transport equation (2.22) there are two types 
of "collisions," one due to Nambu-Goto evolution and another due to interactions. As was 
already stressed, the Nambu-Goto "collisions" take place with the same rate regardless of 
the trajectory and can always be factored out. In contrast, the rate of interactions depends 
on the trajectory in question. Such collisions give rise to an exponential sensitivity to initial 
conditions and contribute to the overall sum of the positive local Lyapunov exponents in 
(6.2). Then the effect of interactions can be estimated from (2.17) as 

/sijB(A,B) apa(A)pfe(B)exp(-7pp | (1, A) A (1, B) A (1, u) A (1, v)| ^"i/^) . (6.3) 

where 7 is some constant and M-i/' is the time scale of collisions. Note that, in general, 
one must solve for the mean directions (i.e. u = (A) and v = (B)) to determine all of the 
parameters of the distribution. However, if we are only interested in a leading-order effect of 
the interactions, then we can assume that the true SRB distribution is not too far from an 
equilibrium distribution described by (5.4) and (5.5). By treating the nonequilibrium effects 
as a small perturbation, we obtain 

/siiB(A,B)aexp^^ + ^-7PP | (1, A) A (1, B) A (1, u) A (1, v)| ^"^/^^ . (6.4) 

This distribution together with equations (4.7), (4.8) and (4.9) can, in principle, be used to 
determine an improved set of hydrodynamic equations. 



7 Summary 

The main objective of the paper was to develop a hydrodynamic approach to cosmic strings 
based on the ideas of the kinetic theory originally introduced in Ref. [21]. In what follows 
we will summarize the key results. 
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First of all, we derived a simplified version of the transport equation (2.22) with Nambu- 
Goto evolution, interactions (such as reconnections of nearby strings and production of string 
loops) and background gravitational effects (such as Friedmann expansion) taken into ac- 
count. Secondly, we proved an i7-theorem (3.4) and derived the four equivalent equilibrium 
conditions for strings (3.16). We also proved a conservation theorem (4.6) and derived a 
system of conservation equations (4.7), (4.8) and (4.9). 

Finally, we derived a complete set of hydrodynamic equations (5.12), (5.13) and (5.14) 
and discussed their possible generalizations in the context of the nonequilibrium steady states 
described by a Sinai-Ruelle-Bowen distribution (6.3). These equations provide an enormous 
simplification to the problem of string dynamics, but we will not attempt to solve them 
here. Instead, we will finish by noting that on a Minkowski background and in the limit of 
small dispersions a and /3, the linearized equations take a particularly simple form (5.18), 
(5.19) and (5.20). In fact, they are not much more complicated than the standard hydro- 
dynamic equations for particles (i.e., continuity equation, heat equation, and Navier-Stokes 
equation), and there should not be any obstacles to solving these equations numerically, if 
not analytically. 
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